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Accurate characterization of charged impurity centers is of importance for the electronic devices 
and materials. The role of valley-spin degeneracy on the screening of an attractive ion by the mobile 
■ carriers is assessed within a range of systems from spin-polarized single-valley to six-valley. The 

screening is treated using the self-consistent local-field correction of Singwi and co-workers known 
as STLS. The bound electron wave function is formulated in the form of an integral equation. 
Friedel oscillations are seen to be influential especially in two-dimensions which cannot be adequately 
accounted for by the hydrogenic variational approaches. Our results show appreciable differences at 
. certain densities with respect to simplified techniques, resulting mainly in the enhancement of the 

impurity binding energies. The calculated Mott constants are provided where available. The main 
conclusion of the paper is the substantial dependence of the charged impurity binding energy on the 
valley-spin degeneracy in the presence of screening. 



C3 



o 
o 



oo 

o 



73.20.Hb, 71.55.-i, 71.45. Gm 



I. INTRODUCTION 



Impurities are introduced to real electronic systems either intentionally as in the form of doping or unintentionally 
mainly in the growth process. Ionization of these impurities contribute charged impurity centers to thej-system and 
^ ■ these_play important role in device operation and material properties; for reviews, see Bassani et ala and Ando 
et al.a for the three-dimensional (3D) and two-dimensional (2D) systems respectively. Due to its fundamental and 
technological importance, the binding of electroaSpto charged impurity centers in the presence of screening has long 
been investigated by the researchers both in 3DtT~H and 2E)|13 _ lL3. Most of these works in 3D aimed to deal with the 
metal-insulator transition usingEJlla Mott's initial approachpwhich he devised, in-fact, for the stability of the metallic 
phase. These approaches can be grouped as variationa]dLnl3~Ej and numericalErEl treatment of the bound electron 
wave function. The former has appealing simplicity but may not be suitable to use in thj^ problem as the resultant 
binding energy will inevitably be higher than .the true ground-state energy. Martino et alu have shown this to be the 
case by comparing Krieger and Nightingale'scl variational hydrogenic wave function treatment with their numerical 
method. Consequently, variational techniques are refined by replacing the trial wave functions with the Hulthen's 
form,E3 resulting in a better agreement with Ref. |?]. In all these works several iprms of dielectric screening have been 
employed such as Thomas- FermiJjTJ Random Phase Approximation (RPA),H~Q and. Hubbard-ShamE'ElS In our work 
we use the so-called self-consistent local-field correctioi*-of Singwi and co-workersEZI to be referred to as STLS. This 
technique is one of the best improvements of the RPA,I13 considering both-the Pauli and Coulomb holes-around the 
r* ■ electrons that take part in the screening of any longitudinal disturbance S3 Very recently Borges et altB also dealt 
with impurity binding energies in 3D using STLS dielectric function, with Hulthen variational wave function. 

The aim of this work is primarily to assess the role of valley and spin degeneracy on the screening and hence binding 
strength of charged impurity centers both in 2D and 3D. The valley and spin degeneracy means, a number of conduction 
band minima and two spin states are energetically degenerate in the absence of symmetry-breaking perturbations like 
strain and magnetic field. The need for this work stems from the current status of electronic devices, evolved into two 
classes as: multi-valley systems dominated by silicon and germanium and single-valley systems realized using mainly 
gallium arsenide. Another source of motivation for the present work-is related to the receatly observed Metal- Insulator 
transition in a 2D system, Si MOSFET, towards zero temperature.E3~E3 We commented^ that this observation could 
be duetto valley phase transition in Si inversion layers, a many-body effect originally anticipated by Bloss, Sham and 
VinterHj For this reason, the effect of valley degeneracy on several physical phenomena needs farther investigation; 
the screening of charged impurity centers is just one of them. We observe that Friedel oscillationsEfliHj associated with 
the screening of an impurity potential are influential in the impurity binding energies. Appreciable differences are 
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seen between our numerical approach based on the computationally-efficient solution of an integral equation and the 
widely-used hydrogenic variational treatments. These results will have implications on the accurate characterization 
of impurities as well as excitonsE3 under the presence of free carrier screening. The screened attractive impurity 
potential has two important aspects: the bound-states and the scattering cross section. In this paper we focus on 
the former and do not address the equally-important problem of the effect of valley-spin degeneracy on the mobility, 
limited by the screened ionized impurities. 

In Sec. II the theoretical details on the variational expressions, the dielectric screening and the integral equation 
formulation are included. The results for 3D and 2D cases can be found in Sec. III. We observe that the phenomena 
in 2D is quite interesting and thereby requires more discussion. Finally, we conclude in Sec. IV. Throughout the text 
we mention the assumptions and simplifications made. 



II. THEORETICAL DETAILS 



We assume that the mobile carriers are over a neutralizing positively charged continuum, that is, the so-called elefe 
tron liquid (EL) model, where electron-electron interactions are rigorously considered using the dielectric formulationEa 
but ionic lattice and disorder effects are ignored. We introduce a degeneracy parameter for the constituent electrons of 
the EL, and thereby consider the spectrum ranging from spin-polarized and single- valley EL to_s.ix- valley EL to assess 
these exchange effects due to valley-spin degeneracy. We consider one singly-ionized impurityO being immersed into 
an EL, and investigate whether this impurity can trap an electron with the screening of the EL present. The binding 
ability in this model is controlled by three effects: i) the attractive bare Coulomb interaction of the ion that enhances 
the binding, ii) the kinetic energy of the bound electron that tries to overcome the binding and iii) the screening 
of the bare interaction by the free carriers of the EL that weakens the binding of the electron. A negative energy 
bound-state may not be possible if the last two effects win over the first as the concentration of the free carriers in the 
EL is increased. We observe that the competition among these effects becomes more interesting in the 2D case. The 
zero intercept density of the binding energy is traditionally known as the Mott value, which renders the comparison 
of different approaches quite simple; we display our results for the Mott constant in Sec. III. Our single impurity 
consideration poses one important restriction on the density of impurity centers; basically, if the bound electron wave 
functions overlap due to large impurity concentration, then the discrete bound-states broaden into bands by means 
of tunnelling between neighbouring sites. However, the aim in high speed electronic devices is to avoid large amount 
of impurities along the mobile carrier paths. 

We work at zero temperature, and for the 2D case, aiming for general results, we assume np_extension along the 
third dimension (i.e., strictly 2D), where electrons still interact with Coulomb 1/R potentiaI.EII We mainly use 3D 
effective Rydbergs (Ry*) for the energies and denote them with an overbar. Also we introduce length-related reduced 
variables by scaling with the effective Bohr radius, a* B and denote them by the subscript r throughout the text such 
as a r = a/a* B , q r = a* B q, for variables having length and reciprocal length dimensions respectively. 



A. Variational Expressions 

For comjjJeteness we first list the expressions for the hydrogenic variational approach. In 3D case Krieger and 
NightingaleB used a variational approach for the bound electron wave function based on a hydrogenic Is trial wave 
function as 
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with a being the variational parameter, whereas Panat and ParanjapeEl used the same form in 2D with radial variable 
r being replaced by the polar radial variable p as 
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The corresponding binding energies are given as 
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in 2D case, where bare interaction is added and subtracted to achieve faster decaying integrand as suggested in Ref. 
|i~0| . In these expressions 630 and €2D represent the static dielectric screening of the EL, which can be computed at 
different levels of sophistication; in the next section we describe the dielectric function we employ. 



B. Dielectric Function For A General Degeneracy Factor 



The ultimately important quantity in screening is the static dielectric function of the EL. Here we would like to 
introduce an extra label (z/),_to the free carriers of the EL in addition to spin (a) and wave vector (fe) labels, to account 
for the extra valley freedomtil. Then, the zerotbu-order polarization insertion diagramed 7T°(g) is modified as in Fig. [I], 
where G° refers to noninteracting propagator .l3 The bare Coulomb interaction being independent of spin and valley 
labels, suggests the introduction of an overall degeneracy factor gd as gd = g s g v where g s and g v are the spin and 
valley degeneracies respectively. For the spin-polarized single-valley EL we have ga — 1 and for the normal-state EL 
having g v energetically degenerate valleys, gd — 2g v . In turn, 7T° is trivially affected by the overall degeneracy factor 
gd as a coefficient in front. The 3D static dielectric function becomes 
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In this expression and throughout the text q n refers to a wave number normalized to Fermi wave number Uf and 
r s = l/a^dnnsD) 1 / 3 , with n 3 D being 3D free electron density. q r and q n are related in 3D as q n = q r r s (2gd/^) ' ■ 
The Pauli and Coulomb holes surrounding the screening electrons are introduced, by the local-field correction G^ LS , 
which needs a self-consistent calculation as described in Ref. [Tt]. Recently, GoldE3 has investigated the effects of valley 
degeneracy on the local-field correction. He noted that the many-body effects are important even at high electron 
densities (i.e., low r s values). We find Gold's work very useful, however, Gold constrained the local-field-correction 
to a Hubbard- like formEj that led to simplicity in the computation. Based on our previous observations ,Lj we avoid 
this simplification and use the standard approach. 
The static dielectric function in 2D is of the form 
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for q n > 2 



r^STLS 
U 2D 



(qn) 



where r s in this case is related to 2D electronic density, n2D by r s — l/a* By /Trn,2D- The relation between q n and q r in 
2D is q n = q r r sy /g~d/2. Again the 2D local-ficliLcorrection, Gf^ L5 (q n ) needs to be self-consistently determined at each 
value [©f r s . We refer to available literatureoE3 for details, however, we would like to mention the work of de Freitas 
et al.^Ii which greatly simplifies the labour in the static structure factor calculation by the so-called Ioriatti-IsiharaES 
transformations. We applied their recipe to a 2D EL having an arbitrary spin and valley degeneracy. 



C. Integral Equation Formulation 

In contrast to simplicity of the variational techniques, they must be used with care in problems such as the 
binding energy, where the variational energy only yields an upper bound for the true ground-state energy. As 
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a better alternative, we present below an approach that leads to an integral equation for which we also develop 
computationally-efficient operator techniques. 



1. Formulation for 3D 



The bound electron feels a centrally symmetric radial potential, and for the lowest Is state, Schrodinger equation 
in 3D becomes 
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UzD,scr{ r r) is the screened potential energy due to a singly-ionized attractive impurity in real space to be computed 
as 
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where r r is the reduced distance in real space; again we add and subtract unscreened Coulomb potential for compu- 
tational reasons.EHl 

The principal problem in the numerical solution is the infinite domain of the wave function. As a remedy, Martino 
et alu noting the difference with the hydrogen atom problem due to the presence of screening, set UsD.scrM to zero 
for distances greater than some large value R. Then, for r > R the wave function for bound-states becomesa 
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upto a normalization constant. The continuity of the wave function together with its derivative at r = R, or 
equivalently, the continuity of the logarithmic derivative of the wave function yields 
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Eq. (Q) is a two-value differential equation problem dealt with shooting type numerical techniques .0 We, rather, 
prefer to convert the radial Scrodinger equation to an integral equation as 



r(r r ) = S(r r ) - dr' T(r') 
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The energy eigenvalue is determined from T(R r ) = — R r y \Eq \ — 1. In our work we extracted the energy eigenvalue Eq 
from the above equation by sampling the wave function in the-5% neighborhood of R r = 10. Eq. ( |ll|) is a nonlinear 
integral equation of the Volterra type in the fixed-point form.c3 However, we observed very slow convergence of the 
standard techniques; for this purpose we first express the Eq. (fTO) as an operator equation as 
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and resort-to operator form of the Newton's method0 which requires the inverse operator of the derivative (Frechet 
derivative-^!) of the operator P evaluated at the function T(r r ). This inverse operator acting on T(r r ) is given as 
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here the prime on the left hand side designates a derivative, whereas, the primes on the right hand side are used to 
produce dummy variables. We retain the first two terms and approximate the final equation as 

lWr r ) = r n (r>) - P [T n (r r )] + — [' dr' T n (r')P [T n {r')] , (15) 
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where T„ denotes T at the n iteration; we determine T„ + i by mixing T new and r„. The final form, then offers a 
rapidly converging algorithm, once we initiate the process at low densities (like r s — 20) using the variational wave 
function as the initial guess and gradually increase the density. 



2. Formulation for 2D 



2D Schrodinger equation for the ground-state wave function reads 

-J + (Eq - U 2 D,scr(Pr)) ^o(Pr) = 



j_ d / ^#0^, . + rE Q -u 2D , scr ( Pr )) MPr) = o. ( Hi) 



p r dp r \ dp 

The expression for the screened potential energy due to a singly-ionized attractive impurity in real space, which shows 
Friedel oscillations is 
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where p r is the reduced 2D radial coordinate and J is the zeroth-order cylindrical Bessel function of the first kind. 
See Fig. || for the screened potential energy at several values of the 2D electronic density; also note the evolution of 
the Friedel oscillations as the density decreases. As in 3D, we work with the function T(p r ) — p r dlnipo/dp r , rather 
than with the wave function itself; in this way an exponentially decaying function is mapped to a linearly decreasing 
one. The nonlinear integral equation satisfied by T becomes 

T( Pr ) = s(p r ) - / dp'^L, (is) 
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where 



S(p r ) = dp'p' U 2DtScr { P ') + ^(^i, (19) 



which is to be computed with very high precision. A nonlinear equation needs to be solved for the bound-state energy 
eigenvalue, of the form 



T(R r ) = -i? rV ^Xi(i? rV ^)/Ko(i?rvW): (20) 

where K is the modified Bessel function of the second kind. 

To achieve much faster convergence than the fixed-point form, the operator P is introduced as 

P [TV)] = T(p r ) - S(p r ) + r dp' 1 ^ = (21) 

Jo P 

We give the final form of the iterative equation we use in 2D which closely resembles the 3D case 

r new ( Pr ) = T n (p r ) - P [T n (p r )} + 2 / " dp'^^P \T n (f/)] . (22) 
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III. RESULTS 



A. 3D EL Results 

We investigate the binding energy of the impurity electron as a function of the electron density and valley-spin 
degeneracy. In Table B, we list the so-called Mott constant as a function of the degeneracy parameter gd from spin- 
polarized electrons to six valley degeneracy as in the conduction band of silicon. Here Mott constant is defined as 

1/3 

a *B n 3D > wnere "3D is the density at which the binding energy reaches zero. It can be seen that for gd greater than 
4 the exchange effects do not lead to appreciable changes in the Mott constant. The spin-polarized EL (gd = 1) has 
the highest Mott constant, which is due to poor screening of the impurity potential by the participating electrons 
having large Pauli holes around them. In Fig. || we plot the variational (hydrogenic) and integral equation solutions 
for the binding energy of the normal-state single-valley (gd = 2) EL; the deviation is clearly visible towards the Mott 
constant. The behaviour for other valley-spin degeneracies is the same apart from a translation according to the Mott 
constant value; refer to Table |. 

Our treatment is based on an isotropic effective mass for the screening electrons, however, mass anisotropy is pre- 
dominantly effective in multi-railcy materials such as silicon and germanium. We refer to available works considering 
the mass anisotropy problem.ErQo In single-valley systems the conduction band effective mass is close to isotropic 
such as the GaAs or the Al^Gai_ x As system. Our calculated Mott constant value for this system (i.e., gd — 2) is 
0.23. Gold and Ghazalicl have very recently dealt with the 3D impurity binding energies using STLS type screening 
and numerical solution for the bound electron wave function. They reported for the same 3D Mott constant the value 
0.25. We attribute the difference between our and their results to the fact that these authors enforced Hubbard-like 
form for the local-field correction which differs from the exact STLS local-field correction leading to a discrepancy in 
the dielectric function. 



B. 2D EL Results 

In general, dimensionality is effective in almost all electronic properties; for our concern in 2D the role of Friedel 
oscillations-|is-enhanced (see Fig. ^|). However, quite commonly the in-plane wave function for (quasi) 2D bound 
impuritiesliH~t2l in the presence of free carriers (i.e., screening) has been chosen to be of e~ p / A type where A is the 
variational parameter. In Fig. ^| the hydrogenic variational probability distribution is compared with that of the 
integral equation solution. The screened attractive potential energy is also added in this figure to aid the comparison. 
The probability distribution obtained by integral equation solution is lower in the first repulsive part of the potential 
energy and higher in the neighboring attractive region than the variational solution; in turn, the electron is expected 
to be more tightly bound. This is seen to be the case in Fig. [| showing the 2D impurity binding energy for the spin- 
polarized and normal-states respectively. Furthermore, a critical Mott density does not exist for these two cases and 
negative energy bound-states are available for all densities, unlike the 3D case. For gd = 4, the variational approach 
predicts a density range, r s — 0.38-1.81, where the binding energy vanishes. The integral equation solution suggests 
that this window is narrower and situated around r s = 1 as can be seen in Fig. ^. As the valley-spin degeneracy further 
increases to gd = 8 and 12, the binding energy curves resemble those in 3D cases; the corresponding Mott densities 
are at r s = 1.52 and 1.48 respectively based on the integral equation solution whereas the variational approach leads 
to higher values (see Fig. ^) . From these three figures we can also conclude that the hydrogenic variational technique 
is successful for the small values of the degeneracy factor, gd- For a better comprehension, we combine all different 
degeneracy cases in Fig. ||, and present a larger density range, till r s = 10. The main observation in 2D is the 
remarkable dependence of the charged impurity binding energies on the valley-spin degeneracy in the presence of 
screening. 

Fig. U shows an interesting strengthening of binding at the high density limit, r s — > 0. The screened interactions 
for several values of gd at r s — 0.02 are plotted in Fig. [|. In this limit, the Friedel oscillations diminish and the 
screening is determined by the exchange effects. Hence, for the spin polarized case (gd — 1), the screening electrons 
cannot approach to the ion due to their sensible Pauli holes resulting in poor screening of the ion potential and 
enhanced binding. As the degeneracy parameter, gd is increased till 12, the influence of the Pauli hole is weakened 
and the central well region of the screened impurity potential gets narrower so that negative energy bound-states are 
no longer supported. Fig. [To] compares the effect of several dielectric functions (RPA, Hubbard, and STLS) on the 
impurity binding energy for gd =2. Sizeable quantitative differences are observed, and STLS dielectric function is 
seen to have stronger screening power leading to a weaker binding. Note the agreement of the three for 7' s — > as 
expected. Hubbard dielectric function follows STLS at the high density end where exchange effects are dominant. 
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Finally, in Fig. || we observe that for gd = 4 case RPA result becomes even qualitatively different and gives negative 
energy bound-states for all densities. 



IV. CONCLUSION 



We investigate the role of valley-spin degeneracy on the screened charged impurity centers. Several complications 
are suppressed for easy comprehension and computational simplicity. These include the mass anisotropy, the effects 
of disorder and ionic lattice on the mobile carriers, the finite temperature, overlap of neighboring bound-state wave 
functions (impurity band formation) and the finite well-width in the 2D case. However, screening is treated using 
the STLS self-consistent local-field correction scheme and the bound electron wave function is handled numerically 
without resorting to simplistic approximations. We observe that care in these two points is rewarding, proven by the 
appreciable differences as compared to widely- used RPA and variational techniques, respectively. We anticipate that 
similar conclusions can be drawn for the Wannier excitons in the presence of free carriers screening. Recently, Ping 
and Jiangc3 investigated the effect of screening on the exciton binding energy in GaAs/ALcGai-xAs quantum wells 
using a rather simple approach based on the Debye screening model and a variational-perturbation method for the 
binding energy. Therefore, the present analysis merits to be extended to excitons. 

The dependence on valley-spin degeneracy is very significant, especially in 2D. From the current electronic devices 
point of view, Si-based and GaAs-based devices are shown to have marked differences in the behaviour of screened 
charged impurity centers. For GaAs-based devices Pauli exclusion principle is more influential in the screening and 
impurity binding energies are larger than in Si-based ones. Binding energy dependence on the degeneracy parameter 
gradually saturates both in 2D and 3D for gd > 8. Finally, the transport through screened charged impurities is also 
expected to have high sensitivity to the valley-spin degeneracy. 
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TABLE I. Critical Mott density, r sc of the 3D EL. The corresponding Mott constants, defined as a B n 3D are indicated in 
parentheses. The numerical results based on the integral equation solution are more reliable (see text). 



Degeneracy factor: gd 


1 

(spin polarized) 


2 

(single-valley) 


4 

(2-valley) 


8 

(4-valley) 


12 

( 6-valley) 


Variational-Hydrogenic: r sc (Mott Const.) 
Numerical: r sc (Mott Const.) 


2.25 (0.275) 
1.44 (0.430) 


3.57 (0.174) 
2.70 (0.230) 


4.11 (0.151) 
3.62 (0.171) 


4.13 (0.150) 
3.80 (0.163) 


4.05 (0.153) 
3.75 (0.166) 



G° a , v (k) 




G o , v (k+q) 



FIG. 1. Zeroth-order polarization insertion diagram, n°(q) for a fermionic system having spin (a), valley (y), and wave 
number (k) labels. 




FIG. 2. Potential energy distribution due to a screened, singly-ionized attractive impurity versus distance. A normal-state, 
single-valley EL is considered (g d — 2) at several densities. 




FIG. 3. Binding energy of a bound impurity electron within a normal-state, single-valley {go, = 2) 3D EL versus the Mott 

1/3 

constant, defined as a B n 3D . Solid line refers to integral equation solution which gives a lower energy than the variational 
treatment based on the hydrogenic wave function denoted by the dashed lines. 




FIG. 4. Probability distribution of the bound electron wave function within a 2D EL having r s = 1 and Qd = 2. Solid line 
is based on the integral equation solution and dashed line refers to 2D hydrogenic wave function, i.e., Eq. (|2j). Also shown by 
dotted lines is the screened potential energy experienced by the bound electron. 



0.0 





FIG. 5. Binding energy of a bound electron within a 2D EL versus r s for <?d=l and 2. Calculations are based on the integral 
equation solution (solid lines) and 2D hydrogenic variational wave function (dashed lines) . 




FIG. 6. Binding energy of a bound electron within a 2D EL versus r s for Pd=4 (i.e., two valley degeneracy). Solid line denotes 
the integral equation solution and the dashed line denotes 2D hydrogenic variational wave function result, both utilizing the 
STLS screening. Dash-dot line refers to RPA screening based on the integral equation solution. 
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FIG. 7. Binding energy of a bound electron within a 2D EL versus r a for <7<j=8 and 12 (i.e., four and six valley degeneracies). 
Solid lines denote the integral equation solutions and the dashed lines denote 2D hydrogenic variational wave function results. 
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FIG. 8. Binding energies in 2D EL versus r s based on the integral equation solutions. For comparison purposes several g d 
values are included. 
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FIG. 9. Potential energy distribution due to a screened, singly-ionized attractive impurity versus distance. The effect of the 
degeneracy factor g d is illustrated from spin-polarized (g d — 1) to six valley degeneracy (gd = 12); all at a very high density 
(r s = 0.02) of a 2D EL. 




FIG. 10. The effect of dielectric function on the binding energy for 2D EL using RPA, Hubbard and STLS screenings; all 
computed by solving the integral equation. Also the STLS screened binding energy is shown based the 2D hydrogenic variational 
wave function labeled by STLS(H). 



